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ABSTRACT 

§3709 S S' 

The electronic oscillator strengths, f, for the N 2 second positive ^ 
+ 

and N 2 first negative systems were determined by measuring the radia- 

* 

tion from the shock layer about a hypervelocity projectile with a time- 
of -flight scanning spectrometer. The measured values at the Av = 0 
bands were f = 0.057 ±0.01 for N 2 ( 2 +) and f = O.O 53 ±0.01 for N 2 (l-). 
The variation of f with wavelength for both systems was found to be 
similar to that previously reported by Nicholls . 




INTRODUCTION 


Measurements of radiant emission from shock -heated gases have 

(1) 

been demonstrated to yield values of oscillator strengths in reason- 

(2) 

able agreement with values obtained by electron beam excitation. In 
previous studies of radiation from shock-heated gases, the gas was 
energized by shock waves generated in shock tubes. The present investi- 
gation utilizes the gas in the shock layer behind the bow shock waves 
of small blunt models flying at hypersonic velocities into nitrogen as 
the source of radiation. Continuously recorded spectra are obtained 
photoelectrically by means of a time -of -flight scanning spectrometer 
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and reduced to obtain oscillator strengths (f-numbers) for the N2 
second positive C 3 it u - B 3 jtg and first negative - X 2 £g 

systems. Nitrogen was chosen as the test gas because of its relatively 
well-known thermodynamic and spectroscopic properties as well as its 
importance in atmospheric physics and as a contributor to the radia- 
tive heating of vehicles entering planetary atmospheres at hypersonic 
speeds . 

EXPERIMENT TECHNIQUE 

The models were lauched from the two-stage shock -heated 0.28- 
caliber light-gas gun of the pilot hypervelocity free -flight facility 

( 3 ) 

at the NASA., Ames Research Center . v The experimental setup is 
sketched in Fig. 1. The models were made of aluminum or a polyfor- 
maldehyde plastic. Because of the short flight distance the aero- 
dynamic heating of the models was such that the aluminum surface remained 
relatively cool, whereas the surface of the plastic model ablated. How- 
ever, as there was no measurable difference in the data obtained with 
the two different model materials, the plastic models were used more 
extensively since they were easier to launch successfully. In addition 
to the time -of -flight scanning spectrometer (described below) the 
instrumentation included shadowgraph stations and electronic counters 
to obtain the velocity and attitude of the model and an Abtronics 2HS 
image -converter camera. Relatively broad -band radiometers, each con- 
sisting of a multiplier phototube fitted with entrance and collimating 
slits and a band-pass filter, were used as a check on the over-all level 
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allowed crude determination of the spatial distribution of radiation 
in the flow field. The use of such broad -band devices is described by 

* ( 4 ) 

Craig and Davy . 

Samples of the nitrogen were taken immediately prior to shooting 
and analyzed with a C.E.C. 21620 mass spectrometer. Typically, the 
contaminants, oxygen, argon, and water vapor, were held to trace amounts 
and never exceeded 0.2 percent. 

TIME -OF -FLIGHT SCAMIHG SPECTROMETER 

Spectra were obtained photoelectrically by means of a time -of - 
flight scanning spectrometer.* This device is shown schematically in 
Fig . 2 . As the model flies by in the focal plane of the collecting 
mirror, the luminous gas cap acts as a moving entrance slit, sweeping 
out the spectrum of the shock -heated gas on the exit slit. (That the 
gas cap is the only source of radiation in the flow field and is indeed 
slit like in form is shown in a typical image -converter photograph of a 
plastic model and in a typical oscillogram of a broad-band radiometer 
(Fig. 3) •) The energy passing through the exit slit of scanning 
spectrometer is divided; 95 percent passes through a splitter plate 
and on to the cathode of an RCA 1P28 multiplier phototube . The output 
of this multiplier phototube is recorded on an oscilloscope and yields 
a continuously recorded spectrum over a wavelength range dictated 
mainly by the geometry of the system. For the tests reported in this 
worl., the wavelength range was as a rule from O.29O to O.L3O p. 

Figure 4 (a) shows a typical oscillogram. The remaining 5 percent of 

*The concept for such a device appears to have originated with E. L. 

(5) 


Woodcock and B. Jones. 
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the energy passes through a narrow-band interference filter (100 A 
half width) and onto a Dumont 6935 multiplier phototube. The output 
of this side band is displayed on an oscilloscope which has been trig- 
gered so as to start sweeping simultaneously with the oscilloscope 
which records the spectra. The side -band oscillogram thus allows a 
specific wavelength to be transposed to the oscillogram displaying the 
spectra. Such wavelength transposition was found to be repeatable from 
shot to shot to ±10A. 

The time -of -flight scanning spectrometer was constructed by 
modifying a Farrand uv-vis f/ 3.5 grating monochromator. A large 
opening was cut in the monochromator housing in the region of the 
entrance slit and the collecting mirror was replaced by a spherical 
mirror which focused at the center of the test chamber. Further, the 
new collecting mirror was masked down to a rectangular slot to improve 
depth of field and to insure uniform illumination of the grating. The 
rise time of the complete system was 25 nanoseconds and was dictated 
by the oscilloscope. 


CALIBRATION 

Wavelength calibration of the spectrometer was accomplished with 

a mercury lamp and appropriate color filters, both on a bench and in 

the test section. The complete optical train was calibrated with a 

tungsten ribbon filament lamp which had been previously compared with 

(g ) 

a primary standard supplied by the National Bureau of Standards . 

The calibration curve for the complete optical train is shown in Fig. 4(b). 
The calibration was checked before and after each series of experiments. 
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EXPERIMENTAL RESULTS 

Time -of -flight scanning spectrometer oscillograms were obtained 
for a range of model velocities of from 4.6 to 5.8 km/sec and ambient 
densities of 0.07 to 0.2 amagat . The calculated equilibrium tempera- 
tures throughout the shock layer ranged from 6000° to 7700° K. A 
typical oscillogram is shown in Fig. 4(a), and Fig. 4(c) shows this 
oscillogram reduced to absolute levels through the calibration curve 
(Fig. 4(b)). The wavelength scale on the oscillogram is nonlinear 
because the model flies on a tangent to the focal plane of the collecting 
mirror and because of slight nonlinearities in the oscilloscope horizontal 
sweep. Results such as shown in Fig. 4(c) were found to be quite repro- 
ducible. The wavelengths of the peaks were repeated to within ±10A over 
the entire range of experiments or at about the limit of accuracy for 
reading the oscillogram. In addition two shots at almost identical 
flight conditions gave results which agreed within 10 percent in absolute 
levels of radiation. 

The peak values of radiant intensity were never greater than 
1 percent of an equivalent black body, calculated at an average shock- 
layer temperature, so self -absorption could be neglected. 

CALCULATIONS 

The electronic oscillator strengths were determined by fitting 
the reduced experimental spectra to theoretical spectra which were 
calculated assuming the electronic oscillator strengths were unity. 

The calculation consists of several steps: first, calculating the 
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radiation per unit volume at every point in the shock layer; then 
integrating over those portions of the shock layer which can be seen 
by the spectrometer. If the radiating gas may be assumed to be both 
optically thin and everywhere in local thermodynamic equilibrium, the 
only 'unknowns are the oscillator strengths of the various radiating 
systems. The details of these steps are described below. 

Radiation per Unit Volume 


There are several method proposed for computing the radiation 

per unit volume from a hot gas a> 9,10 ^ The method chosen for use 

(7) 

here was that of Williams and Treanor, v ' which has the advantage of 
being both straightforward and in a form which allows for extremely 
rapid calculation with a large digital computer. In addition, the 
assumptions and inherent precision of the Williams and Treanor method 
are compatible with the experiments reported herein. 

Williams and Treanor assume that the shape of an electronic - 
vibrational band may be adequately represented by an averaged Q 
branch; that 'is, they choose AT = 0, and take J, the rotational quan- 
tum number, as a continuous variable. This leads to a value of the 
absorption coefficient kyt v " of the form 


k v t v »(v) = C(v , v") exp 


B- 


'V 


hc( V-Vyty") 
kT 


kT J ^ 

where CCv’v") is a quantity, independent of J, whose value is a func- 
tion of the particular vibrational transition v’-v", v is the wave 
number, v° t v » is the band head wave number, and and BJ. are the 


rotational constants for the lower and upper electronic levels, respec- 
tively. The integral of the absorption coefficient over a single line 
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is written 

/ k vV’<iv= g-N 2 f v » v « (2) 

where N 2 is number of molecules in the absorbing state and f y f v " 
is the vibrational oscillator strength. Integrating Eq. (l) gives 

J Vv" dv = C , (v *v") | BJ. - B;i(2J+ l) exp[-B^J(j + l)hc/kT] (3) 

II 2. may be calculated since the density, temperature, and the partition 
functions for rotation, vibration, and electronic excitation are known, 
and thus C(v’v") may be evaluated if fv*v" is known, if Eq. (2) is 
equated with Eq. (3); k v .« v » may then be calculated from Eq. (l) . The 
absorption coefficient at any wave number is then equal to the k v t v " 
summed over all bands which contribute at that wave number . For a band 
system degraded to the violet, 

k = X/ kv ' v " ( k ) 



and the radiation per unit volume from an optically thin gas is 

= k»nB X (5) 

where n is the number density of radiating molecules, is the 
black -body intensity, and k* is the apparent absorption coefficient; 
that is, k T = k[l - exp ( -hcv/kT)] . If more than one band system is 
present, the total radiation at any wavelength is the sum of I-\ for 
each system. 

In the present work the rotational and vibrational constants were 
obtained from Herzberg (ll) and the fraction of molecules in the various 
electronic states was taken from Gilmore (12) . The procedure for 
obtaining fy» v " is discussed in detail below. 
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Radiation Observed 

If the radiating gas is optically thin, the radiation at a given 
wavelength deduced from measurements like those in the present tests 
may be written 

2n 

I^(r,0,cp)F(0,9)r 2 sin 0drd0dcp (6) 

n o o 

Figure 5 shows the coordinate system; F(0,cp) is the fraction of the 
volume not blocked from observation by the model. For a model flying 
with zero angle of attack the flow field is axisymmetric, so 1^ is 
not a function of the azimuth angle cp. Calculating 1^ as a function 
of r and 0 requires a knowledge of the complete density and temperature 
field behind the bow shock. While methods for calculating inviscid 
density and temperature fields are becoming available for certain types 

( 13 ) 

of simple bodies,' ' these methods are, in general, long and expensive, 
and it was felt that such "exact" calculations would add little to the 
accuracy of the final results desired here. Instead the following 
approximate procedure was applied. The shock layer was assumed inviscid. 
The shock envelope was assumed spherical but not concentric with the body 
(cf. Fig. 6), and a correlation suggested by Seiff^ 14 ^ was used to 
calculate the shock standoff distance. Shadowgraphs obtained during 
the tests proved the latter two assumptions quite adequate. A typical 
shadowgraph is shown in Fig. 6. 

If the nitrogen shock tables of Ahtye and Peng ^ 15 ^ are combined 

(is) 

with the tables of species distribution found in Treanor and Logan , v 
the temperature, density, species distribution and, hence, radiation 
per unit volume can be calculated for the gas immediately behind the bow 


E^( watts/p) = 


\R+5(0) pQ 
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shock as a function* of 0 . Further, assuming isentropic flow for the 

dividing streamline and calculating the pressure distributions by the 

(17) 

method of Kaattari permits the radiation from the fluid adjacent 
to the body to be calculated as a function of 0. Then 1^(0) is taken 
as the mean of the values calculated at 0 behind the shock and at 
the body. Since the variation of radiation per unit volume between 
the shock and the body at any given 0 was never more than 20 percent, 
such a linear approximation would seem reasonable. The variation of 
1^ with 0 can be measured by the experimental technique described 
by Givens, et al.^ 18 ^ 

The above calculation procedure was followed at each band head 
until a complete spectrum was obtained over the range of wavelength 
covered by the experiment. 


Oscillator Strengths 

The vibrational or band oscillator strength f v » v " may be written 
in terms of the electronic oscillator strength^ 10 ^ f 


fv‘v" = f 1v'v" 


or the electronic transition momen 


0.(19) 


fv*v" = 


3RoO e S.Q 


where q v , v ii is the square of the overlap integral or Franck-Condon 
factor, R m is the Rydberg constant for infinite mass, e is the charge 
of the electron, and a Q is the radius of the first Bohr orbit. Further, 
f and Rg may be functions of wavelength because of their dependence 
upon internuclear separation. This dependence is usually expressed by 
an experimentally determined variation of R e with the r -centroid (r) 
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/ 'l'v , V' r dr 

r = -k 

r 

J w ^ 

where \|r v i and are the vibrational wave functions of the v* and v" 

levels, respectively, and r is the internuclear separation. 

In order to extract meaningful values of f or R e , it is necessary 
to have accurate values of the Franck -Condon factors. Nicholls^ 20 ^ has 
calculated Franck -Condon factors for the N 2 (2+) and N 2 ( 1-) systems, 
assuming a Morse potential. For the N 2 (2+) system the q v » v » seem 
satisfactory. On the other hand, for the N 2 (l-) system the q v » v » 
calculated using a Morse potential are suspect considering the anoma- 
lous shape of the potential curve of the B state of the N 2 mole- 

cure.*^ 21 ^ For the present work f for the N 2 (l-) system was calcu- 

(22) 

lated using the experimental values of f , v t, reported by Nicholls, 
normalized by his value of f v i v " a "^ bhe 0-0 band and using Nicholls* 
Franck -Condon factor^ 20 ) only at the 0-0 band. This procedure retains 
Nicholls * experimentally determined variation of f with wavelength 
for N 2 (l-) without reliance on the less accurate values of q v i v »- 

RESULTS 

The electronic oscillator strengths at the A v = 0 bands of both 
systems were obtained by fitting the calculated spectra to the experi- 
mental spectra by matching the areas under the radiation per unit wave- 
length versus wavelength curves; that is. 


*1 am indebted to Miss M. Williams for bringing this to my attention. 



where is the absolute experimental intensity. This procedure 

leads to an unambiguous evaluation of the f -numbers at A/ = 0 as the 
Ar = 0 bands of the two systems are not overlapping to any appreciable 
extent. Further, for the N 2 ( 2 +) system, an experimental variation of 
f with wavelength was determined by similarly fitting the remaining 
bands . Figure 7 shows the comparison of the calculated spectra with 
the experimental spectra and it is seen that the' agreement is quite 
good throughout the spectral range covered. 

V ■ A basic assumption in the theoretical calculations used here is 

that the gas in the shock layer is in complete thermo -chemical equi- 

(23) 

librium. It has been found in previous studies both in shock tubes 

• (3) 

and ballistic ranges that under certain circumstances the radiation 
observed may be considerably enhanced by radiation from the nonequilib- 
rium relaxation zone behind the shock front. The possibility of an 
appreciable nonequilibrium contribution was investigated by varying 
the ambient density. In the reduction procedure used here, the non- 
equilibrium radiation would show up as an erroneously large f number, 
the value of which should approach the correct f number as the den- 
sity is increased. Such variation of apparent f number with density 
was discovered and is shown in Fig. 8 , where the apparent f number 
is plotted as a function of the calculated equilibrium density immediately 
behind the bow shock at 0 = 0 . It is seen that the f number of both 
systems levels off to a constant value at densities greater than about 
O .9 amagat, indicating that the assumption of equilibrium is valid 
only at the greater densities. 
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The f numbers at the Av = 0 band determined in this manner 
are f = 0.053 ±0.01 for the N^l-) system and f = 0.057 ±0.01 
for the ^( 2 +) system. The limits of error are estimated from the 
uncertainties in calibration and in the various flow -field approximations . 

Table I compares the f numbers in the present work with those 
reported in the literature. The Bennett and Dalby result has been 
scaled to the 0-0 band as in reference (22) so as to allow direct com- 
parison ■with the present work. The results given here are seen to be 
in quite reasonable agreement with those of other workers despite con- 
siderable variations in experimental technique . This agreement is even 

better than it appears since the shock -tube results have been reported 

(1 9 ) 

as upper limits and any possible systematic error in the electron 

beam experiment due to cascading effects would cause the f number 

reported by Bennett and Dalby to be low. 

The measured variation of the square of the transition moment with 

the r-centroid for N 2 ( 2 +) is shown explicitly in Fig. 9. The agree- 

(22) 

ment with the variation given by Nicholls is good. The points of 

(19) 

the other workers are taken directly from the article of Keck, et al. 

(24) 

and the r -centroids are from Wallace and Nicholls. 

CONCLUDING REMARKS 

The technique of observing the radiating gas cap of a hypersonic 
projectile with a time -of -flight scanning spectrometer has been shown 
to yield quantitative molecular spectra consistent with previous results 
which were obtained with more conventional methods . The advantages of 
' the method used here are very high signal-to -noise ratio, relatively 
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good resolution, and a continuously recorded spectrum.* The main 
disadvantage lies in the lack of complete knowledge of the source of 
radiation, what is, the hypersonic shock layer. As time goes on, 
undoubtedly our knowledge of the details of the flow field within the 
shock layer will improve and thus add to the precision of the method 
reported here. 
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~ *Recently a tec hn ique has been presented using a scanning spectrom- 

( 25 ) 

eter with a shock tube which also yields a continuously recorded spectrum. 
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TABLE I 


Inve st igat ion 

System 

f 

Excitation Technique 

Bennett and Dalby 

N 2 + (l-) 

K 2 (2+) 

0.0438 ± 0.002 
0.0425 ±0.005 

Electron beam 

Allen, et al. 

n 2 + (i-) 

0.09 ±0.05* 

Shock tube 

Keck, et al. 

N 2 (2+) 

0.007 ±0.05 

Shock tube 

Present vork 

r 2 + (i-) 

N 2 (2+) 

0.053 ±0.01 
0.057 ±0.01 

Hypersonic projectile 


*It has been observed^ 20 ^ that a value of f = 0.06 would be a better 
fit to data presented in Ref. (l) . 
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Fig. 9.- Transition probability as a function of 


r-centroid for N 2 (2+). 
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ABSTRACT 


The electronic oscillator strengths, f, for the N 2 second positive 
and N^ first negative systems were determined by measuring the radia- 
tion from the shock layer about a hypervelocity projectile with a time- 
of -flight scanning spectrometer. The measured values at the Av = 0 
bands were f = 0.057 ±0.01 for N 2 (2+) and f = O.O 53 ±0.01 for N^Cl-)- 
The variation of f with wavelength for both systems was found to be 
similar to that previously reported by Nicholls . 

INTRODUCTION 


Measurements of radiant emission from shock -heated gases have 
been demonstrated^ to yield values of oscillator strengths in reason - 

( 2 ) T 

able agreement with values obtained by electron beam excitation. In 
previous studies of radiation from shock-heated gases, the gas was 
energized by shock waves generated in shock tubes. The present investi- 
gation utilizes the gas in the shock layer behind the bow shock waves 
of small blunt models flying at hypersonic velocities into nitrogen as 
the source of radiation. Continuously recorded spectra are obtained 
photoelectrically by means of a time -of -flight scanning spectrometer 
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and reduced to obtain oscillator strengths (f -numbers) for the 
second positive C 3 jt u - B 3 n g and N^ first negative B 2 ^ - X 2 S g 
systems . Nitrogen was chosen as the test gas because of its relatively 
well-known thermodynamic and spectroscopic properties as well as its 
importance in atmospheric physics and as a contributor to the radia- 
tive heating of vehicles entering planetary atmospheres at hypersonic 
speeds . 

EXPERIMENTAL TECHNIQUE 

The models were lauched from the two-stage shock -heated 0.28- 

caliber light -gas gun of the pilot hypervelocity free -flight facility 

( 3 ) 

at the NASA., Ames Research Center . v ' The experimental setup is 
sketched in Fig. 1. The models were made of aluminum or a polyfor- 
maldehyde plastic. Because of the short flight distance the aero- 
dynamic heating of the models was such that the aluminum surface remained' 
relatively cool, whereas the surface of the plastic model ablated. How- 
ever, as there was no measurable difference in the data obtained with 
the two different model materials, the plastic models were used more 
extensively since they were easier to launch successfully. In addition 
to the time -of -flight scanning spectrometer (described below) the 
instrumentation included shadowgraph stations and electronic counters 
to obtain the velocity and attitude of the model and an Abtronics 2HS 
image -converter camera. Relatively broad -band radiometers, each con- 

slits and a band -pass filter, were used as a check on the over -all level 
of radiation seen by the scanning spectrometer. The radiometers also 
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allowed crude determination of the spatial distribution of radiation 

in the flow field. The use of such broad -band devices is described by 

•" (4) 

Craig and Davy . 

Samples of the nitrogen were taken immediately prior to shooting 
and analyzed with a C.E.C. 21620 mass spectrometer. Typically, the 
contaminants, oxygen, argon, and water vapor, were held to trace amounts 
and never exceeded 0.2 percent. 

TIME -OF -FLIGHT SCANNING SPECTROMETER 


Spectra were obtained photoelectrically by means of a time-of- 
f light scanning spectrometer.* This device is shown schematically in 
Fig. 2. As the model flies by in the focal plane of the collecting 
mirror, the luminous gas cap acts as a moving entrance slit, sweeping 
out the spectrum of the shock -heated gas on the exit slit. (That the 
gas cap is the only source of radiation in the flow field and is indeed 
slitlike in form is shown in a typical image -converter photograph of a 
plastic model and in a typical oscillogram of a broad-band radiometer 
(Fig. 3) •) The energy passing through the exit slit of scanning 
spectrometer is divided; 95 percent passes through a splitter plate 
and on to the cathode of an RCA 1P28 multiplier phototube . The output 
of this multiplier phototube is recorded on an oscilloscope and yields 
a continuously recorded spectrum over a wavelength range dictated 
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work, the wavelength range was as a rule from 0.290 to O.43O p. 

Figure 4(a) shows a typical oscillogram. The remaining 5 percent of 

*The concept for such a device appears to have originated with E. L. 

(5) 


Woodcock and B. Jones. 
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the energy passes through a narrow -band interference filter (100 A 
half width) and onto a Dumont 6935 multiplier phototube . The output 
of this side band is displayed on an oscilloscope which has been trig- 
gered so as to start sweeping simultaneously with the oscilloscope 
which records the spectra. The side -band oscillogram thus allows a 
specific wavelength to be transposed to the oscillogram displaying the 
spectra. Such wavelength transposition was found to be repeatable from 
shot to shot to ±10 A. 

The time -of -flight scanning spectrometer was constructed by 
modifying a Farrand uv-vis f/ 3.5 grating monochromator. A large 
opening was cut in the monochromator housing in the region of the 
entrance slit and the collecting mirror was replaced by a spherical 
mirror which focused at the center of the test chamber. Further, the 
new collecting mirror was masked down to a rectangular slot to improve 
depth of field and to insure uniform illumination of the grating. The 
rise time of the complete system was 25 nanoseconds and was dictated 
by the oscilloscope. 


CALIBRATION 

Wavelength calibration of the spectrometer was accomplished with 
a mercury lamp and appropriate color filters, both on a bench and in 
the test section. The complete optical train was calibrated with a 

tungsten ribbon filament lamp which had been previously compared with 

(e) 

a primary standard supplied by the National Bureau of Standards . 

The calibration curve for the complete optical train is shown in Fig. 4(b) . 
The calibration was checked before and after each series of experiments. 



EXPERIMENTAL RESULTS 


Time -of -flight scanning spectrometer oscillograms were obtained 
for a range of model velocities of from 4.6 to 5*8 km/sec and ambient 
densities of 0.07 to 0.2 amagat. The calculated equilibrium tempera- 
tures throughout the shock layer ranged frpm 6000° to 77 OO 0 K. A 
typical oscillogram is shown in Fig. 4(a), and Fig. 4(c) shows this 
oscillogram reduced to absolute levels through the calibration curve 
(Fig. 4(b)). The wavelength scale on the oscillogram is nonlinear 
because the model flies on a tangent to the focal plane of the collecting 
mirror and because of slight nonlinearities in the oscilloscope horizontal 
sweep. Results such as shown in Fig. 4(c) were found to be quite repro- 
ducible. The wavelengths of the peaks were repeated to within ±10A over 
the entire range of experiments or at about the limit of accuracy for 
reading the oscillogram. In addition two shots at almost identical 
flight conditions gave results which agreed within 10 percent in absolute 
levels of radiation. 

The peak values of radiant intensity were never greater than 
1 percent of an equivalent black body, calculated at an average shock - 
layer temperature, so self -absorption could be neglected. 

CALCULATIONS 

The electronic oscillator strengths were determined by fitting 
the reduced experimental spectra to theoretical spectra which were 
calculated assuming the electronic oscillator strengths were unity. 

The calculation consists of several steps: first, calculating the 
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radiation per unit volume at every point in the shock layer; then 
integrating over those portions of the shock layer which can be seen 
by the spectrometer. If the radiating gas may be assumed to be both 
optically thin and everywhere in local thermodynamic equilibrium, the 
only unknowns are the oscillator strengths of the various radiating 
systems. The details of these steps are described below. 

Radiation per Unit Volume 


There are several method proposed for computing the radiation 

per unit volume from a hot gas > ( 7 > 8 -’ 9 - >1 °) The method chosen for use 

(t) 

here was that of Williams and Treanor,' which has the advantage of 
being both straightforward and in a form which allows for extremely 
rapid calculation with a large digital computer. In addition, the 
assumptions and inherent precision of the Williams and Treanor method 
are compatible with the experiments reported herein. 

Williams and Treanor assume that the shape of an electronic - 
vibrational band may be adequately represented by an averaged Q 
branch; that is, they choose AJ = 0, and take J, the rotational quan 
turn number, as a continuous variable. This leads to a value of the 
absorption coefficient ky. t Y « of the form 


k v » v "(v) = C(v'v") exp 


C 


B- 


v 


By-B 


where CKv’v") is a quantity, independent of 
tion of the particular vibrational transition 
number, v° i v " is the band head wave number. 


hc( v-Vy » v ") 
kT 


( 1 ) 


J, whose value is a func 
v’-v", v is the wave 
and and B^ are the 


rotational constants for the lower and upper electronic levels, respec 
tively. The integral of the absorption coefficient over a single line 
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is written 

f I W'^= g-N z f v . v " 

where is number of molecules in the absorbing state and f v » v " 

is the vibrational oscillator strength. Integrating Eq. (l) gives 

J Vv ,,d v =C , (v , v")|B^ - | (2 J + 1 ) exp [-B^J ( J + l)hc/kT] 


( 2 ) 


(3) 


may be calculated since the density, temperature, and the partition 
functions for rotation, vibration, and electronic excitation are known, 
and thus CCv’v") may be evaluated if f v *v" is known, if Eq. ( 2 ) is 
equated with Eq. ( 3 ); k^^" may then be calculated from Eq. (l) . The 
absorption coefficient at any wave number is then equal to the k v : v " 
summed over all bauds which contribute at that wave number. For a band 
system degraded to the violet. 


k = ^ k v v 
v <v 


(4) 


and the radiation per unit volume from an optically thin gas is 

I-\ = k’nB^ (5) 


where n is the number density of radiating molecules, B^ is the 
black-body intensity, and k* is the apparent absorption coefficient; 
that is, k* = k[l - exp ( -hcv/kT)]. If more than one band system is 
present, the total radiation at any wavelength is the sum of for 

each system. 

In the present work the rotational and vibrational constants were 
obtained from Herzberg (ll) and the fraction of molecules in the various 
electronic states was taken from Gilmore (12). The procedure for 
obtaining f v *v" is biscussed in detail below. 
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Radiation Observed 


If the radiating gas is optically thin, the radiation at a given 
wavelength deduced from measurements like those in the present tests 
may be written 

\R+5(q) p0 p 2it 


I^(r,©,cp)F(e,cp)r 2 sin 0drd0dcp (6) 


E^( watts/p) = 

O -rp W V' 

R oo 

Figure 5 shows the coordinate system; F(0,cp) is the fraction of the 

volume not blQcked from observation by the model. For a model flying 

with zero angle of attack the flow field is axisymmetric, so 1^ is 

not a function of the azimuth angle cp. Calculating 1^ as a function 

of r and 0 requires a knowledge of the complete density and temperature 

field behind the bow shock. While methods for calculating inviscid 

density and temperature fields are becoming available for certain types 
(13) 

of simple bodies, ' ' these methods are, in general, long and expensive, 

and it was felt that such "exact” calculations would add little to the 
accuracy of the final results desired here. Instead the following 
approximate procedure was applied. The shock layer was assumed inviscid. 
The shock envelope was assumed spherical but not concentric with the body 
(cf . Fig. 6), and a correlation suggested by Seiff^ 14 ^ was used to 
calculate ' the shock standoff distance. Shadowgraphs obtained during 
the tests proved the latter two assumptions quite adequate. A typical 
shadowgraph is shown in Fig. 6. 

If the nitrogen shock tables of Ahtye and Peng^ 15 ^ are combined 

(is) 

with the tables of species distribution found in Treanor and Logan , 
the temperature, density, species distribution and, hence, radiation 
per unit volume can be calculated for the gas immediately behind the bow 
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shock as a function of G . Further, assuming isentropic flow for the 

dividing streamline and calculating the pressure distributions by the 

. ( i7 ) 

method of Kaattari permits the radiation from the fluid adjacent 
to the body to be calculated as a function of Q. Then I^(e) is taken 
as the mean of the values calculated at Q behind the shock and at 
the body. Since the variation of radiation per unit volume between 
the shock and the body at any given Q was never more than 20 percent, 
such a linear approximation would seem reasonable. The variation of 
1^ with 0 can be measured by the experimental technique described 
by Givens, et al.^ 18 ^ 

The above calculation procedure was followed at each band head 
until a complete spectrum was obtained over the range of wavelength 
covered by the experiment . 


Oscillator Strengths 


The vibrational or band oscillator strength f v » v " may be written 


in terms of the electronic oscillator strength 


(io) 


or the electronic transition moment 


fv*v" — i'lv’v" 
(19) 


R<= 


f v «v" = 


3R0 


i Re 


ea r 


9.v*v 


«ir" 


where q T m is the square of the overlap integral or Franck -Condon 
factor, is the Rydberg constant for infinite mass, e is the charge 
of the electron, and a 0 is the radius of the first Bohr orbit. Further, 
f and Rq may be functions of wavelength because of their dependence 
upon internuclear separation. This dependence is usually expressed by 
an experimentally determined variation of R e with the r -centroid (r) 
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r - J 


f tv»V 


V v j v v "r ar 


t Y ‘V' to 


vhere f v t and tv" are the vibrational wave functions of the v* and v" 
levels, respectively, and r is the internuclear separation. 

In order to extract meaningful values of f or R e , it is necessary 
to have accurate values of the Franck -Condon factors. Nicholls^ 20 ^ has 
calculated Franck -Condon factors for the ^ 2 ( 2 +) and 1-) systems,, 
assuming a Morse potential. For the N 2 (2+) system the q v » v " seem 
satisfactory. On the other hand, for the ^(l-) system the q » m 
calculated using a Morse potential are suspect considering the anoma- 
lous shape of the potential curve of the B state of the N 2 mole- 

cule.*^ 21 ) For the present work f for the 1-) system was calcu- 

(22) 

lated using the experimental values of ? v t v " reported by Nicholls, 
normalized by his value of f v » v " a t the 0-0 band and using Nicholls* 
Franck -Condon factor^ 20 ) only at the 0-0 band. This procedure retains 
Nicholls 1 experimentally determined variation of f with wavelength 
for W 2 (l-) without reliance on the less accurate values of q v t v » • 


RESULTS 


The electronic oscillator strengths at the Av = 0 bands of both 
systems were obtained by fitting the calculated spectra to the experi- 
mental spectra by matching the areas under the radiation per unit wave- 
length versus wavelength curves; that is. 


*1 am indebted to Miss M. Williams for bringing this to my attention. 
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W' 


^ _ J ZV=Q 


A d A 


Av=o 


E A d A 


where is the absolute experimental intensity. This procedure 

leads to an unambiguous evaluation of the f -numbers at Av = 0 as the 
ZV = 0 bands of the two systems are not overlapping to any appreciable 
extent. Further, for the W 2 (2+) system, an experimental variation of 
f with wavelength was determined by similarly fitting the remaining 
bands . Figure 7 shows the comparison of the calculated spectra with 
the experimental spectra and it is seen that the agreement is quite 
good throughout the spectral range covered. 

A basic assumption in the theoretical calculations used here is 

that the gas in the shock layer is in complete thermo -chemical equi- 

(2.3 ) 

librium. It has been found in previous studies both in shock tubes 

(3) 

and ballistic ranges that under certain circumstances the radiation 
observed may be considerably enhanced by radiation from the nonequilib- 
rium relaxation zone behind the shock front. The possibility of an 
appreciable nonequilibrium contribution was investigated by varying 
the ambient density. In the reduction procedure used here, the non- 
equilibrium radiation would show up as an erroneously large f number, 
the value of which should approach the correct f number as the den- 
sity is increased. Such variation of apparent f number with density 
was discovered and is shown in Fig. 8, where the apparent f number 
is plotted as a function of the calculated equilibrium density immediately 
behind the bow shock at 0=0. It is seon that the f number of both 
syst ems levels off to a constant value at densities greater than about 
0.9 amagat, indicating that the assumption of equilibrium is valid 
only at the greater densities . 
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The f numbers at the Av = 0 band determined in this manner 
are f = O.O 53 ±0.01 for the N 2 (l-) system and f = 0.057±0.01 
for the N 2 ( 2 +) system. The limits of error are estimated from the 
uncertainties in calibration and in the various flow-field approximations . 

Table I compares the f numbers in the present work with those 
reported in the literature . The Bennett and Dalby result has been 
scaled to the 0-0 band as in reference ( 22 ) so as to allow direct com- 
parison with the present work. The results given here are seen to be 
in quite reasonable agreement with those of other workers despite con- 
siderable variations in experimental technique . This agreement is even 

better than it appears since the shock-tube results have been reported 
(19) 

as upper limits and any possible systematic .error in the electron 

beam experiment due to cascading effects would cause the f number 

reported by Bennett and Dalby to be low. 

The measured variation of the square of the transition moment with 

the r-centroid for N 2 ( 2 +) is shown explicitly in Fig. 9 . The agree- 

(22) 

ment with the variation given by Nicholls is good. The points of 

( 19 ) 

the other workers are taken directly from the article of Keck, et al. 
and the r -centroids are from Wallace and Nicholls^ 24 ^ 

CONCLUDING REMARKS 

The technique of observing the radiating gas cap of a hypersonic 
projectile with a time -of -flight scanning spectrometer has been shown 
to yield quantitative molecular spectra consistent with previous results 
which were obtained with more conventional methods . The advantages of 
the method used here are very high signal-to-noise ratio, relatively 
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good resolution, and a continuously recorded spectrum.* The main 
disadvantage lies in the lack of complete knowledge of the source of 
radiation, that is, the hypersonic shock layer. As time goes on, 
undoubtedly our knowledge of the details of the flow field within the 
shock layer will improve and thus add cc the precision of the method 
reported here . 
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TABLE I 


Investigation 

System 

f 

Excitation Technique 

Bennett and Dalby 

Na + (1-) 

1^2(2+) 

0.0438 ±0.002 
0.0425 ±0.005 

Electron beam 

Allen, et al. 

k 2 + (i-) 

0.09 ±0.05* 

Shock tube 

Keck, et al. 

N 2 (2+) 

0.C07 ±0.05 

Shock tube 

Present work 

Na + (1-) 

N 2 (2+) 

0.053 ±0.01 
0.057 ±0.01 

Hypersonic projectile 


( 2S ) 

*It has been observed^ ' that a value of f = 0.06 would be a better 
fit to data presented in Ref. (l). 





TEST SECTION 



I* 

S§ 

8 * 

QO 

5 * 


SB 

z_ 

25 

2s« 

*§ 

z“ 

II 


Fig. 2.- Time- of -flight scanning spectrometer. 
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(a) Typical oscillogram from time-of-f light scanning spectrometei . 

(b) Calibration curve of time-of-f light scanning spec^i ometer . 

(c) Reduced experimental spectra. 


Fig. 4.- Steps in determination of experimental spectra. 
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Typical shadowgraph of hypersonic projectile and assumed shock shape. 
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Fig. 9.- Transition probability as a function of r-centroid for N 2 (2+). 
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